Lenses as an Atom-Photon Interface: A Semiclassical Model 
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Strong interaction between the light field and an atom is often achieved with cavities. Recent 
experiments have used a difTerent configuration: a propagating light field is strongly focused using a 
system of lenses, the atom being supposed to sit at the focal position. In reality, this last condition 
holds only up to some approximation; in particular, at any finite temperature, the atom position 
fluctuates. We present a formalism that describes the focalized field and the atom sitting at an 
arbitrary position. As a first application, we show that thermal fluctuations do account for the 
extinction data reported in M. K. Tey et al, Nature Physics 4, 924 (2008). 
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I. INTRODUCTION 

Strong coupling between the light field and a single 
atom is a cornerstone of quantum optics. It has been 
achieved in a large number of laboratories around the 
world in the last decades [J-Q- The interface of choice 
is usually a high-finesse cavity: as well known, the mode 
density and field amplitude is increased when the cavity 
volume is decreased. Recently, some groups have started 
considering a simpler interface: lenses. Focusing is in 
fact the natural way of concentrating light at a desired 
position, here that of the atom. 

This choice is partly motivated by the perspectives in 
quantum information science. Complex information pro- 
cessing tasks will eventually require to scale up from the 
single quantum system to that of a network of systems. 
At each node, information could be processed on atoms 
or ions, whose states can be carefully manipulated and 
efficiently detected; to transmit information between the 
nodes, one would rather encode it in states of light. A 
crucial ingredient of such a network is the coherent op- 
eration of a matter-photon interface at the single-photon 
level 01 • As a first step towards this futuristic vision, 
strong coupling between atoms and propagating light 
should be demonstrated. Experiments conducted on sin- 
gle molecules Q and single atoms [l^, [HI have shown 
promising results. 

A theoretical analysis of such a focusing co nfig uration 
was done previously in [l^, and rectified in How- 
ever, in this latter paper, a significant discrepancy was re- 
ported between some experimentally measured data and 
the corresponding theoretical prediction. As usual, the 
theoretical model does not match the complexity of real- 
ity in many aspects. Here, we show that the discrepancy 
can be essentially accounted for if one includes the fact 
that the atom is at finite temperature and therefore is 
not sitting exactly at the position of the focus, but is 
rather wobbling around it. While this effect is pretty ob- 
vious, its description is not: one must be able to compute 
the electric field felt by the atom at any point close to 
the focus of the system. Due to the strong focusing, the 
paraxial approximation is not valid, so the well known 
closed form solutions for the electric field cannot be used. 

In section |TT] we develop the necessary formalism. We 



start by defining the propagation problem and solving it 
to find an integral equation for the electric field (|II Ap . 
This can later be solved numerically to obtain the field 
anywhere in the system. Further, we model the atomic 
response (jll B[) and the detection process pi C|) . with a 
special attention to the measurement of extinction that 
was reported in [l^ ■ section lllli we apply our 
formalism to quantify the effects of temperature. We 
show that we can fit the experimental data by assuming 
a temperature T « 185 fiK and discuss the meaning of 
this value. 



II. FORMALISM 

A. Electric field 

The electric field after it passes through the single 
mode fiber has, to a good approximation, the Gaussian 
shape 



Eiaser{.P, z ^ Q )=eELe 



(1) 



where we have used cylindrical coordinates {p,(/),z), eis 
the polarization of the field, is the amplitude, and 
represents the width of the beam. 

We model the lens as the infinite plane z = Q. The im- 
pinging electric field at z = 0^ is transformed according 
to 

E{p,z = 0+)=fEip,z^O~). (2) 
The transformation matrix T has the following form: 
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Input laser: ^ Total Electric field: 




FIG. 1: Schematic view of the setup and the corresponding main equations. The atom is sitting somewhere near the focus of 
a confocal arrangement. 



Notice that the factor of , is required by energy 

flux conservation (cf. for instance Pg 391 of [31 )• Now, 
it remains to propagate the electric field in the confocal 
region. We do so with a Green's function approach: the 
electric field is a solution to the following equation, 

V X (V X £;) - k^E = 0. (5) 

Solving for the Dyadic Green's function defined as 

V X (V X t?) - fc't^ = ¥5(f - f ') (6) 

we arrive at 

E{r)^ dS\^k)'—\ -AQ i5(r'), (7) 

where S is the plane where we know fully the electric 
field, R ~ \r — r'\, and wc have used the long wavelength 
approximation kR « 1. Details of the derivation can 
be found in Appendix [K[ 

Although the integrand in this equation is highly oscil- 
latory, we do not make a stationary phase approximation 
as it gives a form similar to the paraxial approximation, 
and this has been shown to give incorrect results [isj . 
We evaluate this integral numerically with usual Gaus- 
sian quadrature llSl. Il6| and the method of Numerical 
Steepest Descent |17| . 

B. Atomic response 

The *^Rb atom, modeled as a two-level system, is ex- 
cited by a weak classical probe beam driving it at the 
resonant frequency uj. The weak beam assumption allows 
us both to ignore incoherent scattering by the atom, and 
to write a fixed phase relation between the dipole and 
laser electric fields. 

The assumption of resonance just means that, with 
respect to the phase of the field at the location of the 
atom, fa, the electric field produced by the atom should 
have a ^ phase lag. 



Combining these with the usual dipole electric field, 
we arrive at, 

4(f) = Mra)\Ct - ^V>^e'^^*-f), (8) 

where Ea and El are the electric field of the atom and 
laser respectively, $l — § represents the phase lag of the 
dipole and e*is assumed to be the polarization of the laser 
field before the lens, and the dipole. 

C. Detectors 

In the experiments [13, [3 , the measured field is col- 
lected by a collection lens, coupled into a single-mode 
fiber, and is subsequently sent to a photon counting de- 
tector. Therefore, only the overlap between the collected 
electric field and the mode of the fiber will be ultimately 
detected. The calibration showed that essentially all the 
input light is detected at the output in the absence of an 
atom. In view of this, we assume that the mode of the 
fiber, g, is proportional to the initial electric field at the 
2 = 0" plane, i.e. 
1 

g (X Elaserip, Z = 0") 

_ 

(X eElc . (9) 
The intensity recorded by the detector becomes 

Idet oc \{g\ Ecou)\'^ 

cx \{ee ^ I E,m)\^ (10) 

In the rest of this article, we will be concerned with the 
extinction, e, of the input field due to the presence of the 
atom defined as, 

Ilaser ^system 

e = -J , (11) 

where Iiaser is the recorded intensity with no atom 
present in the setup, and Isystem is when a single atom 
is present. 
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III. FINITE TEMPERATURE EFFECTS AND 
COMPARISON WITH EXPERIMENT 

In the experiment, the atom is held in an approxi- 
mately harmonic trap at a particular temperature. Thus, 
the two following major effects can be expected: Doppler 
shifts of the incoming laser field, and polarization mis- 
match due to strong focusing by the lenses. 



A. Doppler shifts 

To a good approximation, the potential of the atom in 
the dipole trap is harmonic and can be written. 



H 



I 2 1 2 



(12) 



where = + and cj^ and arc the trapping 
frequencies in the radial and longitudinal direction re- 
spectively. 

We approximate the state of the atom as a thermal 



state 



Patom 



^^n„«) \m){m\^e 



(13) 



here, \m){m\z is the cigenstate of the harmonic oscillator 
in the z direction with eigenvalue = huJzim + ^); 
while Yin is the projector on the eigenspace associated to 
the eigenvalue of = hajp{n + 1), so Tr(n„) = n -I- 1. 

Since the electric field is propagating in the z direction, 
we need only calculate the rms velocity in this direction. 
The square of the momentum is. 



mhuj. 



Tr(atap) + i) = (_ 



2 J (14) 



with ttz = "^^^ Doppler shift to lowest order is 

given by few, where k is the wavenumber, and v is the ve- 
locity. Using an estimate of lOO^K for the temperature, 
we obtain a value of kv « 800kHz. This doppler de- 
tuning is about 0.02 times the 30 MHz linewidth of the 
transition. This justifies us ignoring the Doppler shift 
and considering only the spatial averaging. 



Classical canonical ensemble 



which the simpler classical treatment is justified. There- 
fore, for any function f{r atom, Patom) which depends on 
the position and momentum coordinates of the atom, the 
average value will be given by 



if {"^ atom, Patom)} { f {^atom, Patom) Pat 



^ /e "B-^ f{fatom,Patom)d^pd^r, 



' f{ratom)d^r, (15) 



where we have used the a classical canonical ensemble 

2 

to describe the system and = ^ + V{ratom) is the 
classical Hamiltonian describing the system. 



C. Comparison with experiment 



The extinction defined in ([Tl]) can be reduced to. 



e = 1 



{Ea I El) 



{El I El) 



(16) 



where Ea and El are the electric fields of the atom and 
laser respectively, on the plane of detection. Since the 
propagation of the field is unitary, and we assume infinite 
detectors, we can move the collection plane to that of the 
second lens in the confocal arrangement, z = 2f. This 
means that, 



eixa, 



= 1 - 



3 ~ 

1 {'^a, ^a) {-^a , '^a) 



(17) 



is the extinction of the probe beam in the atom were at 
{xa-,0, Za), where we have taken the atom to be in the 
x-z plane, and I{xa, Za) and K{xa, Za) are shown in Ap- 
pendix |B] Notice that due to cylindrical symmetry of the 
system, it is enough to consider the atom at {xa,Q, Za). 
Finally, the measured extinction must be the averaged 
extinction, weighted in the canonical ensemble as in (|15p . 
So we have to compute. 



Before we continue the averaging over the atomic spa- 
tial profile, it is worthwhile to reconsider in detail how 
the experiment is done. We know that the atom is ini- 
tially loaded into the dipole trap from a MOT, and MOT 
clouds are typically about 0.1 x T^oppier, which in our 
case is about 20^K. 

This typical temperature is justification enough to use 
a fully quantum model to describe the experiment. How- 
ever, rough estimates of the heating due to finite data col- 
lection time put the temperature in the 100/xK range, in 



{(^{ratom)) ^ J ^ ''^^ <^{ratom)d^r 

= - / pdpdze ''^^^ e(p,z).(18) 

We first plot extinction vs temperature for a particular 
value of the focusing parameter, m = where Wl rep- 
resents the width of the initial gaussian beam, and / is 
the focal length of the lens used. 
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FIG. 2: The above is a plot of extinction vs temperature 
for a focusing parameter u = 0.29, ujp = (27r)56 ± 4 kHz and 
bJz = (27r)7± 0.25 kHz as used in the experiments. The dotted 
hues are the effects of the experimental uncertainties in the 
trap parameters. This gives a prediction of the temperature 
of the atom to be about 185 ± 20/iK. 



Now, following reference [13| , we plot curves of extinc- 
tion vs focusing strength, including the effects of a finite 
temperature of the atom, for 10, 50 and ISbfiK. 
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FIG. 3: (Color online) Plot of extinction vs focusing. The 
curve representing the theoretical maximum for an atom sta- 
tionary at the focus (i.e. T — 0), as well as the experimental 
points (filled circles), are the same as in [l^. The other curves 
are plotted for for T =10, 50 and 185fjK respectively. The 
dotted lines come from the experimental uncertainty in the 
trapping frequencies. 

Figure[3]shows that the experimental point seem to fol- 
low the correct trend if we just include a higher tempera- 
ture. This shows that the temperature must be the main 
cause of the discrepancy between theory and experiment. 
Now we plot for a wider range of focusing strengths to 
show the expected trend for a larger focusing. 
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FIG. 4: (Color online) Plot of extinction vs focusing: same as 
in Fig. [S] but for a larger range of focusing parameters. 



D. The meaning of the temperature 

The previous section shows that, from the simulations, 
we should be expecting temperatures of around 180/iK, 
which is higher than the Doppler cooling limit for *^Rb. 
This requires further investigation, since we are expecting 
cooler atoms. We now estimate the initial temperature 
of the trap to evaluate the quality of the trap. 

A simple estimation for the scattering rate would be to 
compare the difference in the number of counts with and 
without an atom in the focus. This simple estimation 
gives approximately 2500s~^ as the scattering rate. Us- 
ing a data collection time of 140ms, and recoil frequency 
of about 3.8kHz gives the energy increase as, 

AE w 2(2500s"^)(140?«s)/i(3.8kHz) (19) 

where the factor of 2 in the energy comes from consid- 
erations of both the absorption and emission processes. 
Since the atom is in a harmonic potential, with a thermal 
density of the form ()13p . we require, 

AE = hupAup + hw^^An^ . (20) 

The above equation can be solved numerically to give a 
temperature increase of, AT « 42^K. The temperature 
obtained from figure [2] is the averaged temperature over 
the probe time. Thus, assuming a linear increase of tem- 
perature, this gives an estimate of the initial temperature 
to be Tinitiai ~ 160/iK. 

Although this temperature is above the Doppler cool- 
ing limit for ®''Rb, it should be noted that the atom is 
loaded from a magneto-optical trap into the dipole trap, 
which is ~ ImK deep. During the loading of the atom 
into the trap, it is unavoidable that the atom heats up. 
Since this trap depth is an order of magnitude larger than 
the Doppler Cooling Limit, the predicted initial temper- 
ature is well within the expected temperature range. 
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We are thus confident that finite temperature effects of 
the atom accounts for the discrepancy reported in [l3| . 
However, these predictions still have to be verified by a 
suitable temperature measurement, one possibility being 
that described in [l^, due to it's high similarity with 
the setup used. Another possibility would be to perform 
Raman sideband cooling, and inferring the temperature 
from there. The latter would seem the better choice, 
since the atom would be cooled in the process. 



Taking the divergence of (|A2|) . and substituting the re- 
sultant equation back, we obtain, 

-(V^+fc^)^ = (^ + ^)^(r-r') (A3) 

This shows that if we have a Green's function which sat- 
isfies the scalar Hclmholtz equation, ie. 



(V2 + fc2)5(r,r") = (5(?^-r'), 



(A4) 



IV. CONCLUSION 



the Green's function as defined in (|A2|) can be con- 
structed as, 



In this work, we present a further step in the char- 
acterization of the results of [l^, We show that 
the inclusion of finite atomic temperature although in- 
tuitively simple, is not straightforward to implement, 
and numerical simulations seem like the only feasible 
option. After we include this in the description of this 
trap, we compare theoretical predictions and experimen- 
tal data, and show that the difference in the data can 
be attributed to the atomic temperature. After numer- 
ical studies, we obtain an effective atomic temperature 
of about Tinitiai ~ 160/xK before probing of the atom . 
Although this temperature seems high, the discussion of 
section UlI Dl gives us some insight as to the problem, and 
we are probably not even justified to make a comparison 
with the recoil or Doppler temperatures as found in the 
literature, due to a slightly differing definition of tem- 
perature used. Thus, we are unable to draw a definite 
temperature for the atom just from this study, and more 
work to measure and perhaps to cool the temperature is 
required. 
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Appendix A: Electric field 

We set out here to give a more formal derivation for 
the Green's function used in equation ([7]) . This approach 
follows closely (loj . 

For completeness, we repeat the following definitions. 
The electric field at any point is a solution of. 



V X (V X i?) - k^E = 0. 



(Al) 



This equation can be solved with the help of the Dyadic 
Green's function defined as. 



V X (V X t?) - = ¥^(7^ 



(A2) 



(A5) 



We use the Green's function representing outgoing spher- 
ical waves. 



giR) 



4:TTR 



(A6) 



where i? = |r — To use the Green's function, we 
approach in a similar fashion to diffraction theory and 
consider the following difference, 

E' ■ (V' X V' X -fc2)t? - [(V' X V' X -e)E'] ■ ^, (A7) 

where E' = E{r"), and V' denotes differentiation with 
respect to the prime coordinates, and in the second term, 
we do not differentiate ^ . Notice that we have preserved 
the order of the multipliction. Then, using Green's sec- 
ond vector identity, this becomes. 



V' 



(A8) 



where tijk is the Levi-Givita tensor, and we assume the 
Einstein summation convention. However, using equa- 
tions (|Aip and (|A2[) . equation (|A7p is nothing but, 
E'S{f— r '). We thus have, 



E'^S{r-r') = V: 



tijkGjXUrak^'lE' 



£ijkEjeimk^lGm\ 



(A9) 



Since we only know the electric field on a plane, we need 
that the Dyadic Green's function be longitudinal on the 
plane, ie. 



n X 



^(f,f') 



= 0. 



f'es 



(AlO) 



In our problem, the plane S, is the z = plane, and so 
n = el- Noticing that ^{f,f') is symmetric in both its 
arguments and that function resulting from operations 
on the right are still solutions of (jA2[) . Then, defining. 



U 



1 -2elel, 



(All) 
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we can see that, 

f ') = f') - t?(f, Ur')U (A12) 



is a solution to ()A2p satisfying the longitudinal condition 
(jAlOp . Taking the curl of (|Al2p . substituting into (jAO]) 
and integrating, we obtain, 



£;(f) = / dS"(ifc) 
Js 



^ 0|^(r'), (A13) 



R R 



which is equation d?]). 



Appendix B: The form of the extinction 

The extinction is given in equation (|T6l) to be 

2 



e = 1 



1 



{Ea I El) 



{El I El) 



Firstly, 



{El \El) ^ // dp'd<P'p\ELe 



TrElwl 



The next inner product is. 



{Ea I El) 



dS 



{f'Ea}' 



(Bl) where I{xa,Za) and K(xa,Za) are, 
I 



(B2) 



-^^r{Xa,Za)k{Xa,Za) (B3) 



J(a;,z) 



r dp' ^£^{cos 9 p, + l)e'^e-^^^/^^ / 
Jo V cos 6p' Jo 



2^~R2- 



(B4) 



^^(•■Z;a,2a) 



-e L 



CQs9pi{p' — Xa cos(/)') + ixsiiicj)') — i^iC sin(?p' — cos6p'{p' — Xa COS (J)') + ixa sin0')e^*'^ j (B5) 



with 



R = ^Jxl + p'^- + zl-2xaP cos (j)', (B6) 
C = 2f-za and (B7) 
r = ^/xl + C + P'^ - 2p'Ccos2 (j)'. (B8) 



Thus we have, 



e{xa,Za) = 1 - |1 2-r(a;a,Za)A'(a;a,^a)P- (B9) 
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